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Abstract 

We point out that the electric dipole moment of the neutron in the Standard Model 
is generated already at tree level to the second order in the weak interactions due to 
bound-state effects, without short-distance Penguin loops. The related contribution has a 
regular nonvanishing chiral limit and does not depend on the mass splitting between s and 
d quarks. We estimate it to be roughly 10~^^ e-cm and expect a more accurate evaluation 
in the future. We comment on the connection between dn and the direct CP- violation in 
D decays. 



Electric Dipole Moments (EDMs) of elementary particles and of composite objects 
allow one to probe CP and T violation at the fundamental level. The Standard Model 
(SM) of particle physics has room for only two sources of CP violation. The first one roots 
in the strong interaction proper, described by the so called '^-term in the QCD Lagrangian 
m. Described by a dimensionless coupling {}, this term is flavor-diagonal and induces a 
large electric dipole moment of the neutron dn ^'&y<3-10~^^e-cm, which by far exceeds the 
experimental limit unless 'd is extremely small, the notorious strong-CP problem. It is 
then natural to assume that "i? is zero or nearly vanishes due to a symmetry or a dynamic 
mechanism [2j. 

The second source of CP violation lies in the electroweak sector of the SM and is 
expressed through the irreducible phase in the CKM matrix Vjk describing the flavor- 
nondiagonal weak transition amplitudes between quarks: 

= ^JlJ' with = ^ V,k UjTf^dk , r^ = 7^(l-75). (1) 

As a result, it primarily manifests itself as CP violation in decays of kaons, B and, 
potentially, D-mesons. The CKM parameterization of flavor dynamics in the SM has 
been successfully tested in detail over the last ten years in K and B mesons, in particular 
through its CP violating thread. 

Being rooted in the intricacies of flavor- changing transitions, the CKM-CP violation 
leads to extremely small effects in flavor- diagonal amplitudes such as EDMs; therefore 
EDMs in general and the neutron EDM dn in particular are a very sensitive probe of the 
underlying source of CP violation. Yet calculating the observable effect for hadrons in 
terms of the fundamental parameters of the SM is often a difficult task requiring control 
of the complicated hadronic dynamics in the grossly nonperturbative regime. 

Following the success of the quark-gluon picture of the hadronic world and of the 
qualitative understanding of the properties of hadrons based on the constituent quark 
model, the early estimates aimed at calculating the EDMs of quarks as the source of the 
nucleon EDM. It turns out that for quarks the SM prediction is particularly suppressed: 
it emerges first at the three-loop level where an additional loop with at least a gluon must 
be included |3]. On top of this, the quark EDM has to be proportional to the quark 
mass; this yields an additional suppression for the light quarks. The same applies to the 
color dipole moments of quarks considered as the simplest induced CP-odd strong force 
generated through weak interactions at small distances. 

It has been noted a while ago [3] that the severe chiral suppression intrinsic to u and 
d quarks can be vitiated in composite hadronic systems like nucleons. The transition 
dipole moments changing (i-quark into s-quark, electromagnetic or color, are suppressed 
by the strange quark mass m^, and such flavor-changing transition without a quark charge 
change are mediated by the loop-induced short-distance renormalization of the bare weak 
interaction at some level via the so-called Penguin diagrams [S]. A similar mechanism is 
believed to dominate the nucleon EDMs in spite of the suppression associated with the 
heavy quark loops, although it is notoriously difficult to account for the long-range part 
of the strong interactions it depends upon [6]. 
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In this Letter we point out that the SM in fact generates the nucleon EDMs aheady 
without any loop corrections, at the tree level of the second-order weak interaction La- 
grangian. This contribution has no chiral suppression neither vanishes if strange and 
down quarks become nearly degenerate in mass. Being loop-free, the induced EDM does 
not suffer from typically small numerically perturbative factors characteristic to Penguin- 
induced effects in spite of the accompanying logarithmic factors. 

The price to pay is that the effect is inversely proportional to the square of the charm 
quark mass. Since the latter is not excessively large in the hadronic mass scale, the related 
suppression comes out relatively mild and this effect may well dominate dn in the SM. It 
is described by multiquark effective operators involving simultaneously u, d and s quark 



The numeric estimate of dn mediated in this way is presently rather uncertain depend- 
ing on the corresponding nucleon matrix elements. We do not even speculate on its sign. 
Applying the natural counting rules we expect \dn\ between 10~^^ and 10~^''e-cm. At the 
same time we expect that definite estimates can be elaborated in the future owing to the 
special form the effective CP-odd operators assume in the CKM case. 

Since the contribution we focus on does not require loop effects, in the following we 
neglect the well known short-distance renormalization of the weak interaction Lagrangians 
altogether. This simplification makes the reasoning transparent and all the expressions 
compact. Including the gluon- mediated short-distance effects is straightforward, and we 
briefly comment on this. 

Let us note that the similar ideas regarding the role of long-distance effects in dn were 
discussed in the 1980s by I. Khriplovich with collaborators [7j as an extension of the origi- 
nal proposal by Gavela et al. [Ij, and then in the 1990s by I. Khriplovich and A. Vainshtein. 
Unfortunately, the latter analysis was not published and its main message remains largely 
unknown. Reportedly, the authors focused on developing a consistent low-energy effec- 
tive theory for the corresponding CP-odd effects, and therefore the loop renormalization 
occupied an important place in the analysis. We are grateful to A. Vainshtein for the 
information and for useful comments. 

1 Effective low-energy Lagrangian 

Through the neutron EDM we are looking into the CP-odd amplitudes that are flavor- 
diagonal in every quark flavor. We assume that the renormalized effective d-ieim. vanishes 
exactly and the only source of CP nonconservation is the general CKM 3x3 unitary mixing 
matrix Vij describing the interaction of W bosons with quarks. Any first-order term in the 
conventional four-fermion weak interaction Lagrangian with AF = is automatically CP- 
invariant due to its Hermiticity. The observable CP-odd effects appear from the second 
order in Cw and thus are proportional to Gf^, being embodied in 



fields. 




(2) 
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The generalized GIM-CKM mechanism ensures that the CP-odd piece of £2 is finite in 
the local four-fermion approximation. 

Descending to a low normalization point we first integrate out top quark and at the 
second stage, below the bottom quark as well. At the tree level integrating them out 
essentially means that all the terms containing t or 6 fields are crossed out of Cyj. 

Within this setup we have effectively nearly a two-family weak Lagrangian 

Cw = JIJ^ with = Vcs cTf^s + Vcd cTf^d + Vus uTf^s + Vud uT , (3) 

yet not quite so, since the four Vki do not form a unitary matrix; in particular, it is not CP- 
invariant. The phases in the four CKM couplings cannot all be removed simultaneously 
by a redefinition of the four quark fields, as quantified by 

A^Im[/4 = ImV;:KrfV;*,Ks. (4) 

Cu, above contains altogether 16 terms bilinear in V and V* . The product C2 of two 
Cw in Eq. (|2]) already consists of 256 termsQ Most of them change quark fiavors - only 
64 are fiavor- diagonal. Of these most still contain the two pairs of complex conjugate 
CKM factors and therefore are automatically CP-even. Hence a handful of terms are 
relevant to fiavor- diagonal CP: only two different operators are not CP-invariant driven 
by the same CKM product f/4, plus their Hermitian conjugated partners. These non-local 
8-quark operators include both q and q fields for each of the four quark fiavors. This is 
readily understood: the CP-odd invariant A (as well as CP-violation altogether) vanishes 
wherever any single CKM matrix element becomes zero. 

Being interested in the nucleon amplitudes, we need to eventually integrate out the 
charm field as well. Here the distinction between the above two terms becomes important. 
If the two charm fields belong to the same four-fermion vertex in the product Eq. as 
in Fig. [1^, they can be contracted into the short-distance loop yielding, for instance, the 
usual perturbative Penguins. These are the conventional source of the long-distance CP- 
odd effects [U [7] . The loop cannot be formed for the alternative possibility where c and 
c belong to different Fig. [T]d since the charm quark must propagate between the two 
vertices. Such contributions therefore are routinely discarded. 

On the contrary, our interest lies in the latter term: it does not involve short-distance 
loops, and has a single charm propagator, although highly virtual in the hadronic scale. 
Each weak vertex contains a flavorless quark-antiquark pair, but these are light down- 
type quarks d and s and are not contracted via a perturbative loop, instead going into 
the nucleon wavefunction. The corresponding operator is 

^ VcsV:^VudV:, j d^x iT{{dV^c){uV^d\ ■ {cV,s){sV,u)^ + H.c. (5) 

The Hermitian conjugate, apart from complex conjugation of the CKM product, is simply 
the exchange between s and c?, s-H-d. 

^This general reasoning does not depend on the tree level assumption: we can exclude Vtj and Vib 
from the pair products using the CKM unitarity. 
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a b 

Figure 1: Two types of CP-odd terms. Weak vertices must be off-diagonal in flavor, either 
for down-type (a) or up-type (b) quark. Solid dots denote the four-quark vertices. Light lines 
correspond to u, d or s quarks, thicker lines stand for charm. 



As the space separation x in Eq. is of order l/rric, eliminating charm results in a 
local OPE; the expansion parameter fiha.dr/'^c is not too small and we need to keep a few 
first terms. The tree-level OPE is particularly simple here and amounts to the series 



c{0)c{x) 



—5'ix) + -^5\x) lip + \5\x) {ilpf + .. 
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(6) 



valid under the T-product. For purely left-handed weak currents in the SM the odd 
powers of l/rric in Eq. ([6]) are projected out, including the leading l/rric piece. We then 
retain only the 1 / term and obtain the local effective CP-odd Lagrangian 



— I 



2m? 



uds'i 



(7) 



u 



in the last expression the covariant derivative acts only on the s-quark field immediately 
following it. 

To address the electric dipole moments we need to incorporate the electromagnetic 
interaction. One photon source lies in the covariant derivative in the operator Ouds-, which 
includes electromagnetic potential along with the gluon gauge field. It is proportional to 
the up-type quark electric charge -|-|. The corresponding photon vertex is local and is 
given by the Lorentz- vector six-quark operator which we denote as O"^^: 

OZ,, = {uY{l-l,)s) (S7^^7"7.(l-75)t^) {h''{l-l,)u). (8) 
Another, non-local contribution is the T-product of the pure QCD part of Ouds 

Ouds = {uY{l~l^)s) (s7^07.(l-75)c^) {h''{l-l^)u) (9) 
with the light-quark electromagnetic currently The total photon vertex is thus given by 



^In general only the sum of the two terms yields the transverse electromagnetic vertex; however, when 
projected on the dipole moment Lorentz structures they separately conserve current. 
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the effective CP-odd Lagrangian 

(10) 

where is the electromagnetic potential and e is the unit charge. 

In principle, the local and non-local pieces above correspond to distinct physics: one 
has photon emitted from distances of order 1 / while the latter senses charge distribution 
over the 1 / /Xhadr range. The latter usually dominates, however the specifics of the left- 
handed weak interactions in the SM makes them of the same 1 / order. 

An interesting feature of the considered contribution is that it remains finite in the chi- 
ral limit and it does not vanish if d and s quarks become nearly degenerate, at first glance 
contradicting the origin of the KM mechanismH This in fact is fully consistent, since the 
external state, the neutron, is explicitly s d non- symmetric. This highlights the differ- 
ence with the short-distance effects for light quarks where severe GIM-type suppression 
would arise and the EDM is proportional to the powers of the light quark masses. 

The CP-odd operators contain strange quark fields. This means that the induced 
effects would vanish in a valence approximation to nucleon where only d and u quarks 
are active. It is known, however, that even at low normalization point the strange sea 
in nucleon is only moderately suppressed. The large-A'^c perspective on the nucleons 
paralleling the picture of the baryon as a quantized soliton of the pseudogoldstone meson 
field |8] makes this explicit: the weight of the operators with strange quarks in the chiral 
limit is generally determined simply by the operator-specific Clebsh-Gordan coefficients 
of the SU{?)) group. This differs from the perturbative Penguin effects which yield small 
coefficients whenever considered in the truly short- distance regime. The neutron dipole 
moment induced by the operator in Eq. ( ITOj) does not need to vanish even in the quenched 
approximation to QCD. 

2 Matrix elements 

The CP-odd operators 0^^^,— O"]^ and Ouds—0\^g have high dimension which is routinely 
regarded as an evidence for being poorly defined for practical applications. However, these 
particular operators possess very special symmetry properties, including antisymmetry in 
respect to s ^ d, which prohibit mixing with lower-dimension operators, and make them 
a suitable object for the full-fledged nonperturbative analysis. 

The neutron EDM is obtained by evaluating the hadronic operator in Eq. f lTU]) over 
the neutron state. Since is T-odd, the matrix element vanishes for zero momentum 
transfer and the linear in q term describes dn'- 

{n{p+q)\C^_\n{p)) = -d^ q^Un{p+q)a'"''j5Unip). (11) 

•^It would vanish if charm and top become degenerate; considering the cases of degenerate bottom and 
strange quarks, or charm and up makes no sense in this context since it has been assumed as the starting 
point that m^jTOc 3> /ihadr while u, d and s are light quarks. 



0:,s + jd'xiT{0^,M J,"^(x)} -H.c. 
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Neither of the matrix elements are easy to evaluate, although one may hope that just this 
particular contribution may eventually be determined without major ambiguity, including 
the definitive prediction for the overall sign. Clearly only the P-violating part of 0^^]^ 
contributes, but we keep them in the original form for the sake of explicit symmetry and 
compactness. 

The contact operator O"^^ is a product of three left-handed flavor currents; Ouds 
instead of the sd current has a flavor non-diagonal left-handed partner of the quark energy- 
momentum tensor in the chiral limit. Therefore it seems plausible that the required matrix 
elements can be directly calculated within the frameworks like the Skyrme model [9l [8], 
or in its dynamic QCD counterpart [10] derived in the large- A^^^ limit from the instanton 
liquid approximation. This is the subject of the ongoing study. 

In the absence of better substantiated calculations we apply the simple dimensional 
estimates to assess the expected size of (i„. We denote 

{n{p+q)\{uLl'^SL) {sLlal^lpdi) {dLl^UL) - (rf o s)\n{p)) = -2i}Cuds q^u{p+q)a'"''y5u{p) 

(12) 

for the local piece given by the operator O'^^^. The reduced matrix element JCuds has 
dimension of mass to the flfth power. We estimate it as 

where /ihadr is a typical hadronic momentum scale and k stands for the 'strangeness 
suppression' to account for the fact that neutron has no valence strange quarks; k^1/3 
is taken typical guess. 

The estimate for dn depends dramatically on the assumed value of /Xhadr- It is known 
that the typical momentum of quarks in nucleon is around 600 MeV or higher. Yet 
six powers of mass in Eq. (IT^ would come from the product of two local light quark 
currents each intrinsically containing factors N^/ Svr^ when converted into the conventional 
momentum representation. This is illustrated by the magnitude of the vacuum quark 
condensate where such a factor effectively reduces /x^adr down to ~(250MeV)^. 

To account for this essential difference we assign a factor of (0.25 GeV)^ =/i^ to each 
additional quark current in the product, while the remaining dimension will be made of 
the powers of /ihadr- Then this contribution to dn becomes 




where A ~ 3.4-10"^ has been used. A potential enhancement may come from summation 
over the Lorentz indices in the currents. 

The most naive estimates for the dn induced by the non-local piece in Eq. ( fTOj) would 
yield a similar dimensional scaling except that no explicit factor Cc = | appears: the 
dimension of the non-local T-product is the same as of O^^^ itself. Within the more 
careful way advocated above the result is literally different: 

|.„|- . e a|32.,>- « ,2. x . (^)' (5^)'. (15) 
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although numerically is not too far away. 

Alternatively, to account for the specifics of the correlator with the electromagnetic 
current in the non-local contribution we can consider the contribution of the lowest reso- 
nant state, the | nucleon resonance A^(1535) referred to below as N. Denoting 

{n{p')\JU0)\N{p))=pj^Un-f5(T>"'quUj^, {N{p')\Ouds{0)-Ouds{0yHp)) = mKdsUj^Un 

(16) 

we obtain for the Feynman diagrams in Fig. [2] 



(17) 



N 



Figure 2: Nonlocal contribution to dn with the intermediate A^. Solid block denotes the 
CP-odd part of the operator Ouds- 

The operator Ouds has dimension ten, and the similar dimensional estimate reads 

IKidsl ~ 1^ f^lfJ'hi.dr- (18) 

can be estimated from the measured transition N ^n+'j approximating the rate with 
the dipole expression 

1 / 

T{N ^ nj) ^ -a,rnP%M'^ y'J^j ^ (380 ± 180) keV (19) 
yielding (0.34±0.08) GeV~^. Finally this estimate would read 
K|->«eA^..»,..,-^.14.10-e.c„.x«(^)'(5|^).(^ 

This value appears quite consistent with the direct dimensional estimate of the non-local 
contribution, in particular considering the fact that the lowest excited state alone may 
not necessary saturate it. 

Therefore, our estimate for dn in the SM literally centers around 10~^^e-cm although 
even the values 5 to 10 times larger may not be excluded. 
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3 Discussions and conclusions 



So far we have neglected short- distance loop effects. They can be included in the standard 
way. Importantly, the gluon corrections do not change the symmetry properties and 
therefore cannot radically modify the structure of the result. 

The evident impact is renormalization of the individual modifying the strength 
and inducing different color flow. This generates different flavor allocation between the 
brackets in the operators Ouds and O"^^, yet does not change their properties. In par- 
ticular, the overall antisymmetry with respect to s and d is never modified. Qualitative 
difference neither arises from the Penguin-induced operators with right-handed currents. 
Most of these changes affect the conventional contributions belonging to Fig. [T^. 

It turns out the situation neither changes where an external magnetic photon or gluon 
attaches to the QQ loop from a single C^'- the sum of all the diagrams vanishes or is 
proportional to the external light quark field mass. In any case such renormalization 
emerges only at the two-loop level and apparently is too suppressed numerically. 

As a result, the short-distance loops should not radically affect dn obtained already at 
the tree level, unless the latter happens to suffer from accidental numeric cancellations. 
All the potential corrections have counterparts in the gluon-induced renormalization of 
the operators Ouds and O"^^ themselves. Therefore we may a priori expect a moderate 
overall enhancement of dn due to increased size of c_ down from the W^-boson scale; the 
final conclusion would require evaluating the matrix elements accounting, in particular, 
for the alternative color contractions in the operators. 

To summarize, we have described the contribution to the EDM of nucleons propor- 
tional to AG'I'/i^j^^j./m^ that does not require loop effects associated with heavy flavors or 
short-distance gluon corrections. It is finite in the chiral limit for light quarks and would 
not vanish even if and were very close. Free from typically small perturbative loop 
factors it plausibly dominates dn in the SM, yet it is interesting regardless of the precise 
magnitude. 

Within the effective theory describing the SM below the bottom quark scale, charm 
quark must necessarily appear for the CP-violation to show up. Therefore, if we avoid 
perturbative charm Penguin loops, they are replaced by the terms suppressed by powers 
of l/nic. We had encountered similar nonperturbative charm effects dubbed 'Intrinsic 
charm' in the inclusive beauty decays [11]; they were distinct, yet not too significant and 
mattered only as long as high precision was aimed at. 

The situation is different for dn and evidently in an important aspect: here we deal 
with 'tree' nonperturbative charm and it is suppressed by powers of I/itlc vs. powers of 
1 / (2mc) for charm Penguins. Moreover, the absence of the loop factor here apparently 
leads to a notably milder suppression. For instance, no as(2mc) enters as a prefactor. 

The left-handed structure of the charged currents in the SM yields the second power of 
the charm mass suppression and, when considering the photon interaction, the local and 
the non-local contribution of the same order. In this situation the local piece possibly 
dominates; this can be verified having more elaborate estimates of the nucleon matrix 
elements. Although rather uncertain at the moment, we believe the proposed contribution 
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to the dn in the SM can eventually be evaluated without a too significant uncertainty, 
owing to a rather peculiar form of the effective CP-odd operator. 

The chiral properties of the SM forbid the effective CP- violating scalar pion-to-nucleon 
coupling at small momenta, which is often important for the atomic EDMs. This applies 
as well to the terms power-suppressed in l/rric- 

Since any CKM-generated contribution to dn must be proportional to AGjp, the consid- 
ered effect is not too far from the top benchmark AGpfi'^^^j.: it is moderated by I /ml, and 
may be somewhat numerically suppressed through explicitly containing the strange quark 
field, a feature likewise unavoidable for CKM. It also does not generate chiral enhance- 
ment ln/i^3^^j./m^ which would be possible with the right-handed currents. Therefore we 
conclude that dn in the SM is moderately suppressed compared to the intrinsic smallness 
of the flavor- diagonal CP violation built in the CKM ansatz. 

The discussed mechanism generating dn roots in the CP-odd interference of the AC = 1, 
AS' = amplitudes in the usual weak decay Lagrangian. This is precisely what induces the 
mixing-free 'direct' CP asymmetry in D decays. Therefore, New Physics (NP) phenomena 
in the latter would directly affect dn at the stated level [12]. It turns out that, barring 
contrived cancellations, the NP dn enhancement is greater than its effect on the D-decay 
asymmetries. This indicates that, in a certain sense, the particular SM implementation 
of the CKM ansatz nearly minimizes the intrinsic size of dn- 

The increase in the EDM depends on the details of the underlying NP interaction. 
For instance, the amplitudes with the right-handed charm field lead to l/rric scaling of 
dn, and the right-handed light-quark currents tend to generate an overall enhancement. 
The latter also induce the CP-odd nNN coupling. Taking the reported asymmetry in 
D^K^K~ and D— )-7r"'"7r~ at face value |13] we typically find an enhancement in dn by 
a factor of 10 to 100, with the case where the NP amplitude is described by ca^^F'^'^'y^u 
alone leading to dn up to 5 ■ 10~^''e-cm. Such an extreme possibility does not look too 
natural from the perspective of speculating about the underlying flavor theory, though. 
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